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Longitudinal waves and two-stream instability in a 
relativistic plasma 
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We investigate the problem of longitudinal waves in a relativistic plasma, with phase velocity lower than 
the velocity of light. We show that such waves exist at arbitrary particle momentum distribution function 
that falls off rapidly enough at large momenta. Earlier papers dealing with these questions are analyzed in 
connection with the problem of pulsar radiation, and the causes of the unclear and confused views 
advanced in some of these papers are discussed. Relativistic-plasma instabilities due to Cerenkov buildup 
of longitudinal waves by beams of high-energy and low-energy particles are investigated. The quasilinear 
relaxation of a beam of high-energy particles in a relativistic plasma is considered and the main 
regularities of this relaxation are explained. 

PACS numbers: 52.40.Db, 52.35.Py, 52.40.Mj 


1. INTRODUCTION 

Laboratory experiments have stimulated the develop¬ 
ment of a theory of collective processes induced by a 
beam of relativistic electrons in a nonrelativistic plas¬ 
ma. 1 This branch of plasma theory is used also in the 
interpretation of astrophysical phenomena (see, e.g., 
Chap. 3 of Ref. 2). At the same time, interest attaches 
in a number of astrophysical problems to the interac¬ 
tion of a beam of relativistic particles with a relativis¬ 
tic plasma. An example is the problem of radioemis¬ 
sion from pulsars (see Ref. 1 and the literature cited 
therein). 

Since two-stream instability is due to Cerenkov in¬ 
teraction of the particles with the waves, oo=kv (oo,k 
are the frequency and wave number, and v is the par¬ 
ticle velocity), the question of two-stream instability 
in a relativistic plasma is frequently associated with the 
question of longitudinal (potential) waves of such a plas¬ 
ma with a phase velocity lower than that of light, t; ph 
= oo/k< c. Longitudinal waves in a relativistic plasma 
were considered initially by Silin 4 and by him with 
Rukhadze 5 for the case of an isotropic Boltzmann dis¬ 
tribution of the particles in energy. Only waves with 


u )/k>c were observed in the cited studies. Silin and 
Rukhadze have therefore made a general statement that 
waves with oo/k< c cannot exist in a relativistic plasma, 
and consequently Cerenkov interaction of longitudinal 
waves of particles is impossible. The problem of lon¬ 
gitudinal waves of particles is impossible. The problem 
of longitudinal waves in a relativistic plasma was later 
investigated by Tsytovich. 6 He has considered a plasma 
having the same momentum distribution function as Silin 
and Rukhadze, 4,5 but, on the contrary, found waves with 
oo/k < c and calculated their damping decrement due to 
the Cerenkov interaction with particles. This result con¬ 
tradicts the aforementioned conclusion of Silin and 
Rukhadze, but Tsytovich made no note of this contra¬ 
diction and did not explain its causes. 

Interest in longitudinal waves in the relativistic plas¬ 
ma has increased after the discovery of pulsars. The 
question of longitudinal waves with oo/k< c was con¬ 
sidered anew by Tsytovich and Kaplan 7 (without refer¬ 
ence to the earlier paper 6 !), but no longer for an iso¬ 
tropic Boltzmann distribution function, as in Refs. 4-6, 
but for one-dimensional power-law distributions. 
Tsytovich and Kaplan 6,7 have stated that longitudinal 
waves with oo/k < c exist and that Cerenkov generation 
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of such waves is possible. This result, however, was 
critized by Suvorov and Chugunov, 8 who reaffirmed the 
previously noted conclusion of Silinand Rukhadze 4,5 that 
there are no waves with oj/k<c , and reached the conclusion 
that Cerenkov generation of such waves is impossible, i.e., 
that two-stream instability is impossible in a relativistic 
plasma. A compromise point of view is advanced by Ginz¬ 
burg and Zheleznyakov, 3 who noted that the realization 
of two-stream instability in a relativistic plasma de¬ 
pends on the concrete form of the plasma particle mo¬ 
mentum distribution function. That this is a moot ques¬ 
tion is stated also in Sec. 5 of Zheleznyakov’s book, 9 
where it is noted, in addition, that if the phase velocity 
of the plasma waves v ph becomes smaller than the velo¬ 
city of light c, then strong Landau damping is to be ex¬ 
pected on particles with velocities v «c, and it is doubt¬ 
ed whether solutions of the dispersion equation e(co,&) 

= 0 exist at all in this region for an equilibrium plas¬ 
ma (c is the longitudinal permittivity of the plasma). 

As to applications of the theory of instabilities of a 
relativistic plasma to the pulsar emission problem, the 
situation seems here quite ambivalent, with two groups 
of authors adhering to two alternate points of view. 

One group 3,8,9 denies outright or by implication the 
possibility of development of two-beam instability under 
pulsar conditions, and take the keystone to be other 
types of instability. The other group, 10 “ 15 on the con¬ 
trary, base their analysis on the concept of two-stream 
instability. 

It is clear from the foregoing that the question of lon¬ 
gitudinal waves and two-stream instability in a relativ¬ 
istic plasma still remains unclear. We must therefore 
first establish which of the preceding results are cor¬ 
rect and which are in error. Such an analysis is given 
in Sec. 2 of the present paper. In agreement with Tsy- 
tovich and Kaplan 6,7 we shall show, first, that plasma 
(longitudinal) waves with phase velocity lower than that 
of light exist at an arbitrary particle momentum dis¬ 
tribution, 1 second, that these waves, even though t> ph 
< c are self-damping, and third, that they can be excited 
by fast-particle beams. In addition, we shall show that 
two-stream instability can occur also at o o/k<c such 
that the dispersion equation e(o>,&) = 0 for the equilibrium 
plasma has no solutions. 

A quantitative analysis of this group of questions will 
be presented in the following sections. At the same 
time, simple qualitative arguments can be advanced to 
explain the situation. 

It follows from the definition of a relativistic plasma 
that its particles have velocities v close to the speed of 
light c y and its momenta p exceed me {m is the rest 
mass of the particle). The degree of relativism of the 
plasma can be characterized by the dimensionless pa¬ 
rameter p/mc , where p is the characteristic momen¬ 
tum of the particle, so that particles with p»p are 
either totally nonexistent or their number is small 
enough (say, exponentially small). If the wave has a 
phase velocity w/k<c , then the velocity of the par¬ 
ticles that are at Cerenkov resonance with this wave is 
v res = u/k, and their momentum is 

Pres =mC [(<*/w) 2 - 1 ] ' 1/2 (1.1) 


Should waves with u)/k< c be strongly damped because 
of the resonant interaction with the particles? Strong 
damping can be expected at p res ^p. Obviously, how¬ 
ever, this damping should be small enough or be com¬ 
pletely absent if 

Pres>P (1.2) 

because of the small number or absence of particles 
with such momenta. From (1.1) and (1.2) follows, in 
analogy with Tsytovich and Kaplan, 7 an estimate of the 
interval of the phase velocities of weakly damped waves 
with co/fc< c : 

0<1— (o/kc<(mc/p) z . (1.3) 

The analysis in Sec. 2 shows that such waves actually 
exist and that it is precisely this interval which contains 
the weakly damped waves obtained in Ref. 6. In con¬ 
trast, Silin and Rukhadze, 5 as well as Suvorov and Chug¬ 
unov, 8 calculated the permittivity and analyzed the pos¬ 
sible realization of waves with oj/k< c in the opposite 
limiting case 

i—(o/kc>(mc/p) i . (1.4) 

Assuming, in addition, that 1 - u/ck «1, Silin and Ruk¬ 
hadze, 4,5 having demonstrated the absence of waves with 

(mc/p) z <l-(o/ck<i, (1.5) 

have interpreted this result as pertaining to the entire 
interval of the possible values of the difference 1 -u;/c& 
« 1, that it is clear from the foregoing that this is in¬ 
correct. 

The foregoing arguments concerned the possibility of 
accepting natural plasma oscillations by a beam. This 
explains also the interest to the problem of natural plas¬ 
ma oscillations. It must be borne in mind, however, 
that even in the absence of natural oscillations with cer¬ 
tain co /k < c, it does not follow, generally speaking, that 
the beam cannot excite perturbations with such co/k. 

This circumstance is well known in the theory of the in¬ 
stabilities of a nonrelativistic plasma. In particular, 
the instability of a slow monoenergetic beam in a hot 
nonrelativistic plasma is due to the buildup of natural 
oscillations of the plasma (see Ref. 16, Sec. 3.3). 

(These perturbations are natural oscillations of the com¬ 
plete plasma + beam system.) It is natural to assume 
that an analogous instability can take place also in a 
relativistic plasma. This is confirmed by the calcula¬ 
tions that follow. 

In Sec. 2 we present a quantitative analysis of lon¬ 
gitudinal waves with o o/k<c in a relativistic plasma ne¬ 
glecting beam effects. In Sec. 3 we consider Cerenkov 
buildup of such waves by a high-energy beam, while in 
Sec. 4 we consider a similar buildup by a low-energy 
beam. By high-energy beam we mean here a beam in 
which the average particle momentum p 0 is large com¬ 
pared with the average plasma-partide momentum, 
p b »P y while a low-energy beam is such that/> & «/>. In 
either case, both the beam and the plasma are assumed 
relativistic, (p b ,p)»mc. 

In connection with the problem of the low-energy 
beam, we note that according to the condition (1.1) and 
the relation p res ~p h «p such a beam interacts resonant- 
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ly with perturbations whose phase velocity satisfies the 
inequality (1.4). However, as noted above, it is pre¬ 
cisely to this interval of phase velocities that the ex¬ 
pressions for the longitudinal permittivity calculated 
in Refs. 5 and 8 pertain. These results of the cited 
papers are used by us in Sec. 4. 

ITie question of two-stream instability of a relativis¬ 
tic plasma is important for a correct understanding of 
nonlinear plasma processes in such a plasma. On the 
whole, this far-ranging problem is beyond the scope 
of the present article. We, however, will touch upon the 
simplest and at the same time important part of this 
problem, namely that of quasilinear relaxation in two- 
stream instability in a relativistic plasma. This is the 
subject of Sec. 5. The results are discussed in Sec. 6. 

2. LONGITUDINAL WAVES IN A RELATIVISTIC 
PLASMA 

We start from the dispersion equation for the longitud¬ 
inal oscillations 

e(tD, k)~Q (2.1) 

with the following expression for the permittivity 4,5,9 

e(<o, A:)=l--(Op 2 < 7 ~ 3 (o)— kv)~ 2 >. (2.2) 

It is assumed that the perturbations depend on the coor¬ 
dinates and on the time like exp(-ivt + ikz). The direc¬ 
tion 2 is arbitrary in the absence of a magnetic field, 
and coincides with the direction of the intensity vector 
of this field if such a field is present in the plasma. 

The symbol () denotes averaging over the particle mo- 
menta, <(...))= //„( />)(... )dp/n, n = Jf 0 (p)dp is the plas- 
ma density, f 0 (p) is the equilibrium particle momentum 
distribution function (one/dimensional), y- (1 +p 2 / 
m 2 c 2 ) l/2 is the relativistic factor, v = c(l + m 2 c 2 /p 2 )~ l/ 2 
is the particle velocity, m is its mass, u) 2 = 4Tie 2 n/m 
is the square of the plasma frequency, and e is the par¬ 
ticle charge. 

We confine ourselves to a one-dimensional particle 
distribution in momentum (the case of greatest inter¬ 
est for applications to the pulsar problem). We find 
with the aid of (2.1) and (2.2) that the longitudinal os¬ 
cillations have a phase velocity equal to that of light, 
a?A = c at k = k 0 , where k 0 is defined by the relation 
k 2 - 2co 2 (y)/c 2 . In this case a ) = u) 0 = ck 0 . Using this re¬ 
sult, we can obtain by successive approximations (cf. 
Refs. 6 and 7) a solution of (2.1) with wave numbers 
close to k 0 . Assuming +fc 0 +Afc, a> = a> 0 +Au>, where 
Afe and Aw are small quantities, and taking into account 
the result of the zeroth approximation e(u> 0 , k<) = 0, we 
obtain in the first approximation 



The zero subscript in the right/hand side of this equa¬ 
tion means that the corresponding derivatives are taken 
at fe = fc 0 ,co = u> 0 . In this case 2 


From (2.3) and (2.4) follows the dispersion law 

(a=c[k—a(k— ko) ]. (2.5) 


Since a > 0, it is seen from (2.5) that waves with k > k 0 
have a phase velocity lower than that of light c o/k<c 
(at k < k 0 we have w/k > c). By the same token we have 
demonstrated the existence of longitudinal waves with 
a )/k<c at an arbitrary particle momentum distribution. 
We note also that the dispersion law (2.5) agrees qual¬ 
itatively with the one obtained in Ref. 7 for a particular 
form of/ 0 . 

We consider now the limits of applicability of (2.5). 

We recall that we have used the assumption that Afc and 
Aco are small, and recognize that in accordance with 
(2.2) and (2.5) this assumption means that 

k 0 (c—v)^>Ak(c—v)—achk. ( 2 . 6 ) 

Since the dimensionless parameter a is of the same 
order of magnitude as 1 -v/c the condition (2.6) actual¬ 
ly means that Afc«fc 0 . It is also clear that, in order of 
magnitude, the dispersion Eq. (2.5) is valid right up to 

We proceed now to calculate the damping decrement 
of the oscillations with oj/k<c when they interact with 
resonant particles. Taking into account the imaginary 
part of (2.2) and assuming it to be a small quantity, we 
obtain the perturbation decrement r (we assume that 
a) = Rea; +ir): 

r=-Im e/(de/d(x))o. (2.7) 

The denominator of the right-hand side of this equation 
is determined by the first relation of (2.4), while the 
numerator is equal to 


I m e = 


f 6 (co —kv) dp =— 

J dp 


4jtVm / 3 df 0 \ 
k‘ V op k-p r „ ' 


( 2 . 8 ) 


The value of the resonant momentum p Tes is given by 
formula (1.1), and if w(fc) is of the form (2.5) we have 
Pres = ™c( 2a A ft Ao)' l/ 2 . 

From (2.7), (2.4), and (2.8) we get 


r-—_ *° cV L>JL) . (2.9) 

2 n <'y 3 (l+y/c) 3 > V dp I p=p ret 

In particular, in the case of an ultrarelativistic plasma 
with a distribution function in the form/ 0 = (nc/2T) 
exp(-c\p\/T), the expression for T means that (cf. Ref. 6) 


nk 0 c 1 

/ 6fc 0 \ 

96to 2 1 

l A k 1 



( 2 . 10 ) 


where y 0 =T/mc 2 . It is seen that at small Afc/fc 0 the 
damping is exponentially small. 

The general relation (2.7) can be represented also in 
the form of the energy balance of the oscillations (cf. 
Ref. 16, Sec. 2.2): 


(£).-«->(£).• 

('( 1+ t)’)/( 1 '( 1 + t))' (2.4) 


dW , l£\l ! 

-= —co Im e —-— 

dt 


( 2 . 11 ) 


where W k = co(ae/aci>)|£j 2 /87r is the energy of the oscil¬ 
lations, and |eJ 2 87t is the energy of the electrostatic 
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field of the oscillations, hi a nonrelativistic plasma we 
have e = l -a>|/a> 2 , so that o>ae/aw = 2, i.e., the total 
energy of the oscillations is equal to double the energy 
of the electrostatic field. On the other hand, for an ul- 
trarelativistic plasma, y»l, as follows from (2.4), we 
have udc/du) *8(y 2 )/(y), i.e., the oscillation energy is 
much higher than the electrostatic-field energy. 

We discuss now the paper of Suvorov and Chugunov. 8 
Their analysis (see also Sec. 5.10 of Zheleznyakov’s 
book 9 ) is based on the assumption that the investigated 
one-dimensional momentum distribution is the result 
of synchrotron radiation with an initial isotropic distri¬ 
bution in the form F(v ) = Qb(v -c), where W is a nor¬ 
malization factor. This approximation is insufficient, 
since the particle velocity is assumed to be exactly 
equal to that of light. To understand the result of this 
seemingly insignificant difference between v and c, we 
assume the initial isotropic velocity distribution in a 
single-velocity form /= Qd(v - vj , where 1 - vjc «1. 
Just as Suvorov, Chugunov, and Zheleznyakov, 8,9 we 
find that after the transverse energy has been emitted, 
a one-dimensional distribution in the longitudinal mo¬ 
menta p „ is established, in the form 


/0>u) = 


n ( i+p 0 z /m 2 c z )' h 
2p 0 (1 +p ll 2 /mV) ,/l ’ 


0, 


p*£p 0 , 

P>Po , 


( 2 . 12 ) 


where p 0 = mv 0 (l -vl/c 2 )" 1 ' 2 . It is seen from (2.12) that 
the system has no particles with momenta p> p 0 . There¬ 
fore, in analogy with (1.3), there exists a phase-velocity 
interval 

0< l—(dlck<';(mclp 0 ) 2 , (2.13) 


which corresponds to a real value of the longitudinal 
permittivity. On the other hand, Suvorov and Chugun¬ 
ov 8 have assumed p 0 =°°, and in this case the interval 
(2.13) does not exist. 3 


We see thus that the reason why longitudinal waves 
with c o/k < c were not observed in Ref. 8 was the choice 
of an equilibrium distribution function unsuitable for this 
purpose. 

As to the work of Tsytovich and Kaplan 7 they considered 
distribution functions in energy ( E ) in the form 


(x—1 )nE. r 1, E<Ei 
( E+E .)" 10, E>E, 


(2.14) 


where E* and £ max are the average and maximum par¬ 
ticle energy, and k is a dimensionless parameter that 
characterizes the rate of decrease of distribution func¬ 
tion at high energies. In addition, they used 7 a suc¬ 
cessive-approximation method similar to that described 
above, and calculated the damping decrement of the 
longitudinal waves in two cases: at 2< x < 4 and at 
x > 4. At x < 4, however, this method cannot be used, 
so that the corresponding results of Tsytovich and Kap¬ 
lan 7 are in error. This error was in fact the basis for 
the criticism of Suvorov and Chugunov. 8 The latter, 
however, criticized Kaplan and Tsytovich’s results 
pertaining both to x < 4 and x > 4, and this is an error 
on the part of Suvorov and Chugunov. 8 


3. EXCITATION OF LONGITUDINAL WAVES IN 
A RELATIVISTIC PLASMA BY A BEAM OF HIGH- 
ENERGY PARTICLES 


The possibility of buildup of longitudinal waves in a 
relativistic plasma by a beam of high-energy particles 
is clear even from formula (2.9): when (dF/dp) p=Pres 
> 0 the growth rate of the oscillations is positive. The 
corresponding instability is a kinetic variant of two- 
stream instabilities and occurs in the case of a beam 
with sufficiently large momentum spread. On the other 
hand, if the beam-particle momentum spread is not 
too large then, just as in the case of a nonrelativistic 
plasma, development of hydrodynamic two-stream in¬ 
stability is possible. 


To ascertain the conditions for the particular in¬ 
stability, we consider by way of example a Gaussian 
distribution of the beam-particle momenta: 


AW—e-exp 


r_ 0 >-p») 2 i 

L p T b 2 J ‘ 


(3.1) 


Here p b is the average beam-particle momentum, p Tb 
is the “thermal” spread of the momenta in the beam. 

It is assumed that p b >p Tb »mc. It is also implied that 
the total plasma + beam system distribution function con¬ 
sists of two parts: f=f 0 +f b , where f 0 is the plasma dis¬ 
tribution function. The plasma is assumed to be rela¬ 
tivistic, p>mc (p is the average plasma-particle mo¬ 
mentum), and such tha.tp«p b . 

Under these assumptions the dispersion equation of 
the longitudinal oscillations of the plasma + beam sys¬ 
tem takes the form 


l+e o +eb=0. 


(3.2) 


Here e 0 denotes the second term of the right-hand side 
of (2.2), and e b is given by 


G) p{ , 2 =4jie 2 rc 6 /m, y b =p b /mc, yn—PTb/mc, 


(3.3) 




is the probability integral with complex argument, and 
n b is the beam density. 

Equations (3.2) and (3.3) are similar to the standard 
equations of the theory of two-stream instability of a 
nonrelativistic plasma (see Ref. 16, Sec. 3.1). Using 
this analogy, we conclude that the beam can be regarded 
as monoenergetic, and that the instability it excites is 
hydrodynamic if x> 1, i.e., if 


'iTb/'ib<\a-kV b \'( b t /kcttr h '( b VkoC, ( 34 ) 

where T h is the growth rate of the hydrodynamic in¬ 
stability. If the inequality in (3.4) is reversed, then it 
must be assumed thaLthe beam has a large momentum 
spread and that the instability it excites is kinetic. 

To express the condition (3.4) in terms of the equili¬ 
brium parameters of the plasma + beam system we must 
solve the dispersion Eq. (3.2) in the monoenergetic - 
beam approximation. In this approximation it follows 
from (3.3) that 
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e 6 =-(D P 6 2 'V (to-kV b ) 2 . 


(3.5) 


(3.12) 


Assuming that w = a) 0 +Aw,fe = fe 0 +Afe, where a> 0 and k 0 
are defined in Sec. 2, and using (2.4), we reduce the 
dispersion Eq. (3.2) to the form 

Aco—c(l—a) AA:=(Wra 0 ) (koc!2yy b ) 3 ((a—kV b )~ 2 . ( 3.6) 

Here y = ((y 3 )) l/s and the angle brackets (• ••) denote 
averaging over the plasma distribution functions; in 
other words, the quantity y characterizes the average 
relativistic factor of the plasma particles. [At/ 0 
~exp(-c|/>|/T) we have y = 6 l/3 y 0 ,r 0 = T/mc 2 , cf. Sec. 2.] 

Next, in analogy with nonrelativistic-plasma theory 
(see Ref. 16, Sec. 1.5), we introduce the quantity u> de¬ 
fined by the relation 

65—A©—c (1—a) A/e—«— kV b . (3.7) 

From (3.6) we obtain for a cubic equation that yields 
the maximum growth rate of the hydrodynamic instabil¬ 
ity of a monoenergetic beam 



When (3.8) is taken into account, (3.4) means that 

(3.9) 

Yb ' Tlo ’ Y 

This inequality is a generalization, to the case of a re¬ 
lativistic plasma (y»l), of the known condition 
y Tb < (n h /n<) l,2 y b under which a beam in a nonrelativistic 
plasma is monoenergetic. 18 It is seen from (3.9) that 
with increasing relativism of the plasma the condition 
for a monoenergetic beam becomes more stringent. 
According to (3.8), when the relativism increases the 
instability growth rate also decreases. Both effects can 
be regarded as a consequence of the increase, noted in 
Sec. 2, of the vibrational energy of the particles of the 
relativistic plasma relative to the electrostatic ener¬ 
gy, a)aeo/aci)~y 2 »l. 

It is known from the theory of the interaction of a 
relativistic beam with a nonrelativistic plasma 18 that 
even a beam with a large thermal scatter can be re¬ 
garded as monoenergetic if (n b /n 0 ) l/3 y b ~ 1. In the case 
of a relativistic plasma, as seen from (3.9), this in¬ 
equality is replaced by the following: 

(Wrc 0 ) 3 Yb>Yo. (3.10) 

We note also that from the second equality in (3.7) we 
obtain the following estimate for the value of A&/fe 0 that 
corresponds to the maximum of growth rate: 

aWo«(y/^) 2 . (3.11) 

By assumption y/y 6 «l, so that the role of the reson¬ 
ant particles of the plasma at such Afe/fc 0 is negligibly 
small. 

We assume now that a condition inverse to (3.9) is 
satisfied, and consider the kinetic instability. The 
corresponding analysis can be carried out starting 
either from (3.2) and (3.3) or directly from the general 
expression for the kinetic growth rate (2.9). As a re¬ 
sult we obtain for the maximum growth rate of the kine¬ 
tic instability the expression 
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It is seen that the relativism of the plasma influences 
the kinetic instability much more than the hydrodynamic 
instability, r kin ~y~ 5/2 (it is recognized that fe 0 ~y l/2 ). 

4. INSTABILITY OF LOW-ENERGY BEAM 

Unlike in Sec. 3, we shall show here that longitudin¬ 
al oscillations of a relativistic plasma can be excited 
also by a low-energy relativistic beam with mc«p b «p. 
We neglect the thermal scatter of the beam. Its per¬ 
mittivity then takes the form (3.3). Assuming c 0 % to be 
a small parameter, we obtain from (3) by successive 
approximations 

(D==/cV' 6 ±0)p 6 (l + eo(©, k)o= h vr>)~ h - (4.1) 

Taking into account our assumption me « p b «p , we 
conclude that the quantity £ 0 (o o,k) in (4.1) should be cal¬ 
culated for u/k satisfying the condition 

1/y 2 <1-g>/&c<1. (4.2) 

According to the statements made in Secs, land2, this 
is precisely the phase-velocity interval to which the cal¬ 
culations in the cited papers 4,5 pertain (see also below 
concerning Ref. 8). We use the results of these papers 
to demonstrate, without further calculations, the effect 
investigated by us. 

1. Plasma with three-dimensional distribution function 
proposed in Ref. 5. 

This function is 

<«> 

where p is the three-dimensional momentum and p»mc. 
With the aid of formula (13.9) of Ref. 5 we find that with 
/ 0 so defined we have 

e„(o>, ^ (in-!—- 2-in) ; (4.4) 

d 2 = T/4ue 2 n 0 . It follows from (4.1) and (44) that at small 
k the instability is of the hydrodynamic type, and at 
large k it is kinetic. The kinetic instability can be in¬ 
terpreted as a dissipative buildup of negative-energy 
waves (see Ref. 16, Sec. 2.3). Hie maximum growth 
rate is obtained at the junction of the hydrodynamic 
and the kinetic instabilities, namely at 

(4.5) 

In order of magnitude we have then 

r^Ctopb/^') (lnY6 2 )' A - (4.6) 

Taking (4.5) and (4.6) into account, we find that the 
beam can be regarded as monoenergetic if 

*Yxfc/'Yfc < ( n b /n Q ) v Y<f ft (In ft 2 ) (4.7) 

where y 0 = T/mc 2 . At a much larger thermal spread of 
the beam, the plasma + low-energy beam system is 
stable. 

We note that in our approximation of large yjy b the 
growth rate turns out to be independent of the degree 
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of relativism of the plasma (of the parameter y 0 ). This 
is due to the character of the dependence of the optimal 
wave number k on y 0 [see (4.5)]. It is also of interest 
to note that the limiting thermal spread of an unstable 
beam increases with increasing y 0 [see (4.7)], and does 
not decrease as in the case of the high-energy beam 
[cf. (3.9)]. 

2. Plasma with one-dimensional distribution function 
of the type (2.12) 


As noted in Sec. 2, such a one-dimensional distribu¬ 
tion function is established as a result of synchrotron 
radiation of the plasma with an isotropic 5-function dis¬ 
tribution in momentum. If the momentum of the par¬ 
ticles with such a distribution is large compared with 
the momentum of the beam particles, p Q »p b , then it 
can be assumed in the calculation of the permittivity 
that Po~* °°* It was in this approximation 
that e 0 (co,fe) was calculated by Suvorov and Chugunov. 8 
Using Eq. (5) of that paper, we have 


e 0 (<o, k) 0= 




(4.8) 


With the aid of (4.1) and (4.8) we find that the maxi¬ 
mum instability growth rate is reached at k= (3 tt) i/2 u p / 
2c, and in this case it is of the order of r ~v Pb /y b . The 
condition that the thermal scatter of the beam be small 
takes the form 


W(AA;) = 


n b mc z k 0 

4a('if 6 -'if m<n )(AA:) 2 ’ 

0, 


Ak(p b )<Ak<Ak(p m{n ) 
Ak>Ak(p min ), Ak<Ak(p b ) 


Pmin/ttlC, 


(5.4) 


Just as in the case of a nonrelativistic plasma, in the 
case of relaxation the beam loses half its initial ener- 

gy . 18 


Following Ref. 1, we consider also the dynamics of 
the quasilinear relaxation of the beam. The distribu¬ 
tion function of the beam is then represented in the form 

# / fn b /(p b -p.(t )), p.(t)<p<p b 

,{p ’ t)= \ 0, p<p.(t) • < 5 - 5 > 

For p*{t) we get the equation 

dX/dt —/tX7(l-X); 

x _ P-(t) n n. WQ>.+0) (5 - 6 > 

p b 2 n 0 Ak b c ’ n W (/>.—0) 


is a quantity of the order of the Coulomb logarithm. 1 
From (5.6) it follows that 


Mt)=pbl[l+{2Aty h ). 


(5.7) 


From this we get the time at which the width of the plat¬ 
eau reaches half its final value (p* = p b / 2): 


1 n 0 Ak 0 c 

t'h*= -j • 

ji n b co P 


(5.8) 


lTbhb< (rib/n 0 ) ' lu ib- (4.9) 

This does not contain the relativistic factor of the plas¬ 
ma, a natural result, since the distribution (2.12) cor¬ 
responds to a moderately relativistic plasma (y> «1, 
although with a long ultrarelativistic tail. 

5. QUASILINEAR RELAXATION OF HIGH-ENERGY 
BEAM IN AN ULTRARELATIVISTIC PLASMA 

We consider now the kinetic stage of quasilinear re¬ 
laxation of a high-energy beam. By the standard method 
(see Refs. 1, 18, 19) we obtain the initial system of 
quasilinear equations 


For a plasma with a momentum distribution function in 
the form/-exp(-c|p| /T) we have a = isy 2 o, k <fi = u P Vo 2 , 
where y 0 = T/mc 2 . In this case 

(5.9) 

3i n b ^6 <jJp 

We see that the relaxation time increases with increas¬ 
ing plasma relativism like yl /2 . 

We can consider in similar fashion spatial rather than 
the temporal relaxation of the beam. 1 However, since 
the group velocity of the oscillations is close to the light 
velocity, it suffices for this purpose to make the sub¬ 
stitution a /a cd/dz. Therefore, 


11 

dt 


5 lJLi f i * 

ac dp 


11 

dp ’ 


(5.1) 


dW 

dt 


fCo \ dp/ P =* Pi h) 


(5.2) 


The oscillation energy W is connected with the electric 
field by formula (2.11), and the relation between the 
resonant wave number and the resonant momentum is 
\k(p) = kj2y 2 . 

In analogy with Ref. 1, we get from (5.1) and (5.2) an 
equation for the quasilinear integral 


dr d k 0 Wm z c i 

rr - 7 — =0 - 

dt L op ap J 


(5.3) 


From this we find that the quasilinear relaxation causes 
the beam distribution function to acquire a plateau in 
the interval Pmn^P^Pb> where p ^ is the scale of the 
average plasma particle momentum. The distribution 
of oscillation energy over the wave numbers takes in 
this case the form 


bu — Ct'i, (5.10) 

is the distance over which the beam loses half its ener¬ 

gy. 


6. DISCUSSION OF RESULTS 

We have demonstrated that in a relativistic plasma 
there can exist longitudinal waves with phase velocities 
smaller than the phase velocity of light, and have es¬ 
tablished that such waves can be exicted via Cerenkov 
resonance by beams in which the particle energy is 
large compared with the plasma-par tide energy. In 
addition, in accordance with our analysis, the long¬ 
itudinal waves can be excited by beams of low-energy 
particles. By the same token we have shown that the 
same variants of two-stream instability exist in a rela¬ 
tivistic plasma as in a nonrelativistic plasma. 

Besides the linear theory of the two-stream instabil¬ 
ities of a relativistic plasma, we have analyzed the 
quasilinear relaxation that accompanies the develop- 
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ment of such instabilities, and obtained quantitative 
characteristics of this process. 

It follows from our analysis that the use of the two- 
stream instability concepts to interpret the radioemis¬ 
sion from pulsars does not contradict the assumption 
that the near-pulsar plasma is relativistic. Investiga¬ 
tions in this direction were performed by a number of 
workers. 10 " 15 It must be noted, however, that in these 
studies the relativism was taken into account by using 
rather crude and not always consistent simplifying as¬ 
sumptions. This raised objections (see Ref. 3). The 
approach described above makes it possible to improve 
the earlier analysis. The development of a consistent 
theory of collective processes in a relativistic plasma 
would be useful both for the interpretation of radiation 
from pulsars and other astrophysical problems, and 
for laboratory applications. 

We are grateful to R. Z. Sagdeev for a discussion of 
the results and to G. I. Melikidze for useful discussions. 


** It is assumed, however, that the distribution function de¬ 
creases rapidly enough at momenta much larger than the 
average; for details see below. 

2) In the derivation of (2.4) we take into account the one-dimen¬ 
sional character of the particle momentum distribution. For 
this reason, expressions (2.4) cannot be obtained from the 
analogous formulas (38) and (39) of Ref. 6, which pertain to 
an isotropic distribution. 

3) The need for introducing a finite p Q is clear also from the 
Maxwellian-distribution example considered in Ref. 17. 
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Nonlinear effects in the excitation of a magnetic field in 
a strongly ionized plasma 

B. A. Al’terkop 

Institute of High Temperatures, USSR Academy of Sciences 
(Submitted 27 July 1978) 

Zh. Eksp. Teor. Fiz. 76, 971-975 (March 1979) 

We consider the nonlinear dynamics of the magnetic instability of a one-dimensionally inhomogeneous 
plasma. It is shown that the principal nonlinear effect that limits the growth of the magnetic perturbations 
is the dissipation of the temperature inhomogeneity of the plasma. An estimate of the maximum value of 
the exciting magnetic field is obtained. 

PACS numbers: 52.35.Py, 52.35.Mw 


1. It is known that to obtain a high temperature dense 
plasma with the aid of powerful photon beams or charged 
particles, the optimal procedure is to produce adia- 
batically spherically symmetrical heating. Therefore 
serious misgivings were expressed when spontaneous 
excitation of a magnetic field, reaching 10 4 -10 6 G, was 
observed in experiments with laser plasma. 1 " 5 The 


appearance of strong magnetic fields can lead to asym¬ 
metry of the heat flux and affect most adversely the re¬ 
gime of compression and heating of the plasma. 

Several assumptions were made concerning the cause 
of excitation of the magnetic field. It is shown in Refs. 
6-8 that the magnetic field can be produced in the re- 
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